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Abstract: We study a Peccei-Quinn invariant extension of the next-to-minimal supersym- 
metric Standard Model (NMSSM), which turns out to be free from the tadpole and domain 
wall problems. Having a non-renormalizable coupling to the axion superfield, the SM sin- 
glet added to the Higgs sector can naturally generate an effective Higgs fi term around the 
weak scale. In the model, the lightest neutralino is dominated by the singlino, which gets 
a mass only through mixing with the neutral Higgsinos. We explore the phenomenological 
consequences resulting from the existence of such a relatively light neutralino. The cou- 
pling of the SM singlet to the Higgs doublets is constrained by the experimental bound on 
the invisible Z-boson decay width. Under this constraint, we examine the properties of the 
SM-like Higgs boson paying attention to its mass and decays. We also demonstrate a UV 
completion of the model in SU(5) grand unified theory with a missing-partner mechanism. 
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1. Introduction 

The next-to-minimal supersymmetric Standard Model (NMSSM) introduces a SM singlet 
S to explain the origin of a supersymmetric Higgs /i term of the MSSM However, if 
5 is a true singlet under all symmetries, it becomes difficult to embed the NMSSM into a 
more fundamental theory such as a grand unified theory (GUT). This is because the GUT 
partners of the MSSM Higgs doublets also couple to S and radiatively generate a large 
tadpole for 5", destabilizing the gauge hierarchy Non-renormalizable interactions are 
another source of large tadpoles A symmetry under which S transforms non-trivially 
can solve the tadpole problem, but generally introduces another problem. If one considers 
a discrete symmetry, dangerous domain walls would be formed in the early universe |5| . 
On the other hand, a global symmetry spontaneously broken by S and the Higgs doublets 
would give rise to an unacceptable visible axion. Additional structure is thus needed to 
generate // dynamically from the coupling of S to the Higgs doublets while providing a 
viable framework for the grand unification. 

In this paper, we point out that the difficulties arising due to the singlet S can be 
avoided in a Peccei-Quinn invariant extension of the NMSSM (PQ-NMSSM) where the PQ 
symmetry is spontaneously broken at a scale much higher than the weak scale by the axion 
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superfield. The PQ symmetry forbids the generation of large tadpoles for S while solving 
the strong CP problem Q. Furthermore, the domain wall problem can be resolved in the 
presence of PQ messengers that couple to the axion superfield. It also turns out that a 
non-renormalizable coupling of S to the axion superfield naturally leads S to get a vacuum 
expectation value around the weak scale. 

The Higgs and neutralino sectors are considerably modified by the addition of S. In 
the PQ-NMSSM, the lightest neutralino consists mostly of the singlino because it acquires 
a small mass through mixing with the neutral Higgsinos after the electroweak symmetry 
breaking. The presence of such a relatively light neutralino leads to phenomenological 
consequences different from other NMSSM models. In particular, the LEP bound on the 
invisible Z-boson decay width places a stringent constrain on the coupling of S to the Higgs 
doublets if the decay mode is kinematically allowed. This constraint becomes important at 
large tan /3. We also note that loops involving the Yukawa coupling of the singlino give an 
additional positive contribution to the SM-like Higgs mass. This contribution is insensitive 
to tan /3, and arises when the Higgsinos are lighter than other MSSM sparticles. Another 
interesting feature is that the decay of the SM-like Higgs boson into a pair of the lightest 
neutralino can be dominant at low tan f3. 

Since interactions of S are controlled by the PQ symmetry, it is possible to embed the 
PQ-NMSSM into GUT models without the tadpole problem. The associated UV comple- 
tion is then related to the doublet-triplet splitting problem. We find that incorporating 
the PQ symmetry in a missing-partner model for supersymmetric SU(5) GUT ||^, ^ can 
lead to the PQ-NMSSM at low energy scales. In addition, it naturally achieves a phe- 
nomenologically acceptable value of the axion decay constant as Fa ~ y/MsusyMpi with 
MgusY being the SUSY breaking scale. The PQ symmetry is also important for suppressing 
harmful dimension 5 operators for proton decays. 

This paper is organized as follows. In the next section, we present the model and 
discuss its general properties and various constraints on the singlet couplings. Then, in 
section 3, we construct a low energy effective theory below the SUSY breaking scale to 
examine the phenomenological aspects resulting from the existence of a singlino-like light 
neutralino. Section 4 is for the discussion on how to UV complete the model. The PQ- 
NMSSM can arise as a low energy theory of a missing-partner GUT model. Section 5 is 
the conclusion. 



2. PQ-invariant extension of the NMSSM 

In this section, we extend the NMSSM^ to incorporate the PQ symmetry and study the 
properties of the model. The PQ-invariant extension turns out not only to provide a 
solution to the strong CP problem [|14| but also to solve both the tadpole and domain wall 
problems. We also examine constraints on the coupling of S to the Higgs doublets. 

^Extensions of the MSSM with a SM singlet have been received revived attention j^, |l^, |ll|, |l^, |l^ after 
the first results on the Higgs search at the LHC had been announced. 
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2.1 Model 



In the PQ-NMSSM, an effective Higgs fi term is generated by the vacuum expectation 
value of S which couples to the Higgs doublets and to the axion superfield X through the 
PQ-invariant interactions 



//" -X"* 
(feXSHuHd+ d^eK——S + h.c., (2.1) 
J Mpi 

for {X, S, Hu, Hj) carrying the PQ charges as (1, 2, —1, —1). All other terms involving S are 
forbidden by the PQ symmetry in the renormalizable superpotential. A superpotential term 
X'^HuHd/Mpi, which is allowed by C/(1)pq, can be removed by a holomorphic redefinition 
of S without loss of generality. Here we have taken such a field basis. 

In the following, we assume that a mechanism to stabilize X is operative with its vac- 
uum expectation value fixed at 10^'^^^^ GeV as required by cosmological and astrophysical 
observations. As suppressed by Mpi, the interactions between X and other scalar fields 
rarely affect the saxion potential at IS"! <^ \X\ <^ Mpi. The Higgs potential can thus be 
examined by replacing X with its vacuum expectation value. Then, the involved mass pa- 
rameters are determined by the SUSY breaking scale MsusY and the axion decay constant 
Fa ~ \X\. We note that, when treating X as a spurion field, the present model can be 
regarded as the nMSSM,^ where the superpotential contains an effective tadpole for S: 

Wes = XSHuHd + ml{l + e^B,)S, (2.2) 

with ~ ^SUSY and 

ml ~ kMsusytt^. (2.3) 

Mpi 

The appearance of ml and terms can be understood by promoting k to a function 
depending on SUSY breaking fields in a hidden sector. It is obvious that the model does 
not suffer from the tadpole problem because a tadpole for S requires a higher dimensional 
couphng of 5 to X as dictated by U{1)pq. For Fa = lO^""^^ GeV, MsusY ~ 1 TeV and 
K less than order unity, the value of rn-o can naturally be around the weak scale. Hence, it 
is natural to expect that electroweak symmetry breaking would occur at the correct scale. 



In fact, the same spirit is shared with the Kim-Nilles mechanism |21] that explains the 
smallness of /Lt in extensions of the MSSM with C/(1)pq. Since S and Hy^ ^ carry ?7(1)pq 
charges and develop vacuum expectation values, the Higgs and neutralino sectors have 
small mixing with X suppressed by Fa- 



^The nMSSM |l6j assumes specific discrete R symmetries to ensure the absence of large tadpoles 
for S. A general discussion on the phenomenological aspects of the nMSSM for small tan j3 can be found 
in p^ , p^ , p^ , p^ . Some cosmological issues have also been discussed in p^ . Neglecting small mixing 
with the axion superfield, the Higgs and neutralino sectors of the PQ-NMSSM have the same phenomeno- 
logical properties as the nMSSM. However, the cosmological properties can be different depending on the 
cosmological evolution of the saxion. A continuous symmetry to restrict couplings of S in the NMSSM has 
been introduced in [EOl, where it is explicitly broken only by a linear superpotential of S. 
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Let us now examine the vacuum structure of the low energy effective theory given 
by ( |2.2D . Including soft SUSY breaking terms, the scalar potential of the extended Higgs 
sector reads 

V=\{9' + 9''){W-\H,\^? + \g'\HiH,\^ 
+ \\HuHd + ml\^ + |A|2|5|2(|i7„|2 + |^^|2) 

+ mlj\Hu\^ + ml^\Hd\'' + + [A^XSH^Hd - B^mlS + h.c.) , (2.4) 

where mf is a soft scalar mass squared, and A\ is a soft ^d-parameter. The potential 
contains four complex parameters. A, Ax, rn^ and -B^- Among them, A and ttiq can be made 
real and positive by a field redefinition of H^^d and X. Furthermore, if sig{A\) = sxg{B^), 
one can rotate away the phases of and i?^ by redefining S. We will assume this is the 
case, for which CP invariance is preserved in the Higgs sector and there is no mixing between 
scalar and pseudo-scalar fields. From the above scalar potential, it is straightforward to 
get the conditions for electroweak symmetry breaking. Similarly as in the MSSM, two of 
them can be written 



mir — 171% tan^ B „ 
2^^" tan2/3-l ^'^^^ 



sin 2/3 = — 7^ 7y — — K 7ir-^, (2.5) 

for /ieff and b^g defined by 

Meff = Xvs, bes = K^xvs + "Iq); (2-6) 

where = vs'mf] and {\H^\) = vcos/3 with v = 174 GeV. The value of |5| at the 

vacuum is fixed as 

AxXv'^sm2/3 + 2B^ml ,„ 
= 2(m| + A2.2) • (2.7) 

The tree- level mass matrices for the scalar fields are presented in the appendix |^. 

To explore the global structure of the potential, one can substitute S by the solution 
of dsV = 0. Then, the Higgs potential (2^) is written 

V-V\ \Ax\HuHd-B,ml\^ 

which increases monotonically along the D-flat direction \H^\ = \H^\ when 

Ri>l and 3i?i > 2 + i?2, or 1 > -R? > R2, (2.9) 

where Ri^2 are defined by i?im| = |(2//eff— ^A)"T'| + (2/Uefr— ^Asin2/3)A2ti2|2/3 and i?2W'5 = 
(2//eff — ^a)^ — + "^l/d ~^ '^t^eS ~ 2^eff)- If the above condition is not satisfied, the 



- 4 - 



potential may develop another minimum away from the weak scale. The involved soft 
parameters are then constrained by the requirement that the electroweak vacuum should 
be a global minimum. For ~ -^SUSY ^eff' which is the case we shall focus on, the 
stability condition ( |2.9| ) requires Aj^ < m'g or m| < Aj^ < m\ with niA being the mass of 
the CP-odd neutral Higgs boson. Keeping this in mind, we will consider also the case with 
Ax ~ MgusY) which is favored to avoid large mixing of the SM-like Higgs scalar with the 
singlet scalar when /Ueg and tan /3 are large. 

An important consequence of f7(l)pQ is the appearance of a relatively light neutralino 
with a large singlino component. This is because the PQ symmetry prevents the singlino 
S from having a supersymmetric mass. The lightest neutralino is mostly singlino if the 
masses of the bino B and wino W are larger than Xv, as is the case for A < 1. The 
neutralino mass matrix for [B, W^, H^, H^, S) is given by 

(Ms -Mz sin 9w cos (3 Mz sin 9w sin /3 \ 

Myi^ Mz cos Ow cos P ~Mz cosOw sinP 

~Mz sin 9w cos /3 Mz cos 6w cos (3 — /^off' ~ Aw sin (3 

Mz sin 9w sin (3 —Mz cos 9^ sin (3 —f^cS —Xv cos (3 

\ -XvsinP -Aucos/3 / 

where 9w is the weak mixing angle. If we write the lightest neutralino as a linear combi- 
nation of (B,W^,Hl^,H^,S): 



,(2.10) 



-° - NfB + N^W° + ivf ° + iVf ° H^ + NfS, (2.11) 



Xi 



then we find 



ivf = i-^(i + ---)4' 

^1° = -eH(sin/3 + • • • ), N^' = -e^(cos/3 + • • • ), 
N? = -eBeH{l + ---), ivf = 6^6^(1 + •••), (2.12) 
for <C 1 and | | ^ 1 . Here the epsilon parameters are defined by 



Xv _ g'v cos 2/3 _ gv cos 2/3 



and the ellipsis indicates terms of higher orders in or can see that ^ 

is needed to make x? massive through mixing. In the following discussion, we will neglect 
small gaugino components of Xi since it does not change our results substantially. Then, 
one can find 

m^o ~ 2 ^^cffiVf "A^f ° + XvN^'^Nf cos /3 + XvN^'^Nf sin /3 

AVs.n2^^^_AV^^/AVXy ^^^^^^ 



^Actually a minimum of the potential does not lie in the D-flat direction unless = wlf^. However, 
as will be shown in the appendix a minimum other than the electroweak vacuum, if exists, is located near 
the _D-flat direction for much of the parameter space. This justifies our approach of examining the D-flat 
direction to see when there can be another minimum. See also [^, where the stability of the electroweak 
vacuum has been examined within the framework of the effective Lagrangian beyond the MSSM. 
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As we will discuss later, the singlino-like neutralino with a small mass can considerably 
change the phenomenological properties of the model. 

Meanwhile, there can exist PQ messengers ^' + ^' which are vector-like under the SM 
gauge group and obtain heavy masses from the coupling in the superpotential. Such 

interaction can play an important role in the saxion stabilization because it induces a 
radiative potential for the saxion after SUSY breaking. The presence of PQ messengers 
also helps to avoid the domain wall problem. Let us consider A^^ pairs of + ^' forming 
5 + 5 representation under SU(5), for which the gauge coupling unification is preserved. 
Then, the domain wall number is given by 



iVow = |Avt -6|. (2.15) 

This implies that the domain wall problem can be resolved for A^^ = 5, 7. For other cases 
with Adw 7^ 1) the formation of dangerous domain walls can still be avoided if the saxion 
is displaced far from the origin after the inflation ends so that the PQ symmetry is not 
restored at high temperatures ||2^ . 



2.2 Constraints on the model parameters 

Since S modifies the Higgs and neutralino sectors, it is of importance to explore constraints 
on the singlet couplings A, Ax, and m|. Here we focus on the case with 0.1 < A < 1 at 
the weak scale as would be natural because an effective fi term is generated as ficS = 
with S fixed around MsusY • Let us first examine the mixing of the singlet scalar with the 
Higgs doublets. After taking the rotation of {H^,H^) by an angle /3, the mass matrix for 
the CP-even scalar fields has 

(M^)ii = M| cos^ 2/3 + X^v^ sin^ 2/3, 
(Ml) 13 = At;(2//efr-^Asin2/3), 

(Ml)33 = m| + AV, (2.16) 

where {Mfj)ii constitutes an upper bound on the mass of the lightest CP-even Higgs boson 
at the tree-level. In the following, we would like to consider the situation that the SM-like 
Higgs boson has negligible contamination from the singlet scalar. For 0.1 < A < 1, this is 
achieved when 



[4s 



^ sin 2/3 



<. 1, (2.17) 



with m|. being of the order of -Mgugy- It is thus found that, if Ax is larger than 2/iefT, the 
mixing can get a sizable suppression at some region of tan /3. One would otherwise need 
/igfj <^ m-l to suppress the mixing. 

In the PQ-NMSSM, a stringent constraint on A comes from the experimental bound on 
the Z-boson invisible decay rate because the PQ symmetry makes Xi light. The singlino 
mixes with neutral Higgsinos to induce the interaction Xi'^'^Xi^n |24], through which Z 
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can invisibly decay into pairs of the lightest neutralino. The coupling for this interaction 
is given by 

9 (^^M\2 ^^THla\ 9 A2i>2cos2/3 



y^XiXi 2cos^H/ V ) ^cosOw /^eff 

where the last approximation is valid for small Xv/fieS- Hence, at large tan/3, the interac- 
tion gets strong while the mass of Xi becomes small. The above coupling mediates the Z 
decay into Xi with 

4y0y? , o / A \ Vcos2/3\ V200GeV\^ 

if Mz > 2m^o. Here = (1 — Am?-Q/M'^)^/'^ is the velocity of x? in the rest frame of Z. 
The process Z — )• XiXi contributes to the invisible Z decay and is tightly constrained by 
the LEP data to occur with a small rate, T^^-^o-j^o < 2 MeV |25]. This translates into 

The above constraint on A around the weak scale becomes important for large values of 
tan/3. To kinematically forbid the mode Z — t- XiXi; we need 

1/2 



which is possible for A < 1 at low tan/3. For instance, at small tan/3 around 2, the LEP 
limits on the Z invisible width exclude values of //eff in the range between 196 GeV and 
260 GeV for A k, 0.6 [|l6|. 

On the other hand, for the theory to remain perturbative up to Mgut, A should be 
small enough at the weak scale. NMSSM models with XSHuH^ usually require A less than 
0.7 — 0.8 for tan/3 > 2. In models with a superpotential term S^, the upper bound on 
A decreases as the coupling for increases. However, the situation in the PQ-NMSSM 
is different because is absent and the PQ messengers with mass Mq, cc \X\ affect the 
running of gauge couplings. Above the scale Mq,, gauge couplings have larger values than 
in the MSSM and slow down the running of Yukawa couplings. This results in the increase 
of the perturbativity bound on A by (5A < 0.1 [Q. A large number of light PQ messengers 
are favored by raising the bound, but disfavored by the requirement of the perturbativity 
of gauge couplings up to Mqut- If there exists an extra gauge interaction, the perturbation 
theory would be valid below Mqut for a larger value of A at the weak scale. 

Though we do not discuss it here in detail, there is also a constraint placed by cos- 
mology. If Xi is the lightest sparticle, its relic abundance should not exceed the measured 
amount of the dark matter. In the present model, the production of dark matter relies on 
the cosmological evolution of the saxion, which has a very flat potential generated after 
SUSY breaking and thus can play some non-trivial role in cosmology. On the other hand, 
the gravitino or axino can be lighter than Xi depending on the mediation mechanism of 
SUSY breaking and on how the saxion is stabilized. 
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3. Low energy Higgs sector 



In this section, we study the low energy Higgs sector. To see the impact of the PQ-NMSSM 
specific Higgs properties, we consider the decoupUng hmit of the MSSM where ah heavy 
Higgs states decouple below MsusY and thus one combination of H^^d behaves exactly like 
the SM Higgs scalar H. In such a situation, we include the singlet S and construct a low 
energy effective theory below MgusY to examine how much the model departs from the 
MSSM. The modification is mainly due to {i) the extra contribution to the Higgs quartic 
coupling, which is a general property of NMSSM models, and (ii) the presence of a light 
neutralino that is singlino-like, which is a consequence of the PQ symmetry. 

3.1 Eflfective theory below the SUSY breaking scale 

For 0.1 < A < 1 at the weak scale, <^ m| is favored to suppress the mixing between 
H and the singlet scalar. Here we consider such a case and assume that the MSSM 
sparticles other than Higgsinos obtain masses of the order of MsusY- The singlet scalar 
is also assumed to have m| ~ M^^^y. For /igflf less than MsusY; the low energy effective 
theory below MsusY contains Hu,d and S in addition to the ordinary SM particles. The 
Lagrangian relevant to our analysis is given by^ 

- C,s = ^{\H\^ - v^f + {yttRQiH'^ + ^^osHuHd + y'uHHuS + y'dH'^HdS + h.c), (3.1) 

where = —ia2H*, and yt is the top- Yukawa coupling. The singlino Yukawa couplings 
at the SUSY breaking scale are 

y'uiMsvsY) = Acos/3, y'diMsvsv) = Asin/3, (3.2) 
while the Higgs quartic coupling is given by 

Ah(MsuSy) = ^ Y + T ^'^ + ^Ai^ltree + ^A^lloop, (3.3) 

where 5\H\tree is the threshold correction coming from tree- level exchange of the singlet 
scalar, and SXh\ioop is from the loops involving the stops: 

AX I ~ A^(2/ieff-^Asin2/3)^ 

dA^/ tree TT^ ' (3-4) 

Zrrig 

'^A^.lioop ^^(^t-^), (3.5) 



^To obtain the couplings of the SM-hke Higgs boson more precisely, one needs to know the mixing 
between the SM-like Higgs boson and singlet scalar. To this end, one can replace the scalar part of ( |3.l[ ) by 

r I Ah — (JAfl-ltrcc „|2 2n2 , ^2 , ,2 2^0 |2 
— l.Gff I scalar = (|-H| -V ) + (mg + A V )\b—Vs\ 

X{2^ca - Axsm2l3) 2 2wo n 

^ ^ ^(l-f^l -V ){S -Vs) +C.C. 

and fieS by AS. The small mixing with the singlet scalar reduces the couplings of the Higgs boson h. The 
reduced couplings can be obtained by taking the replacement 

Au|2^off - ^Asin2^ 

1 ~ o 
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where Xt = {At — /igfr cot /?)^/AfsusY with At being the ^-parameter for HytRQi- 

The physical mass of the CP-even neutral Higgs boson h can be obtained using the 
relation m| = 2Xhv^. For this, we need Xh renormalized at the weak scale. In the effective 
theory, a low energy value of Xh is determined by the renormalization group (RG) running 
equation: 

= li? + 4(3y2 + y'^ + y'^ - ^A)Xn +3B- I2yf - A{y'^ + y'^f) , (3.6) 

with the parameters A and B defined by 

4A = 3g^ + g'^ AB = 3g^ + 2g^ g'^ + g'\ (3.7) 

Here one should note that the mixing between the neutral Higgs boson and singlet scalar 
would slightly modify the running equations. 

To see the qualitative properties of the Higgs mass, we make an approximation taking 
into account that the dominant effects on the RG running come from the term yf, and also 
from the terms y'^^ if A is not small. The Higgs boson mass is found to be approximately 
given by 



for rrit = ytv and fi^s ^ ^h- The first line is the well-known result for the Higgs boson 



mass in the MSSM 28|. On the other hand, those in the second line correspond to the 
additional contributions arising due to S, i.e. as a consequence of the extra Higgs quartic 
coupling A^I-ffu-f^dP; the mixing between the singlet scalar and neutral Higgs boson, and 
the singlino Yukawa interactions affecting the running of the Higgs quartic coupling at low 
energy scales. The last two contributions are approximately estimated as 

2 



A negative contribution to nih from the mixing with the singlet scalar is present in any 
NMSSM model. In the PQ-NMSSM, <C m| leads to a large suppression of this effect 
for Ax < /ieff- For Ax > 2fief[, small mixing can still be obtained at some values of tan/3. 
One should also note that there is a PQ-NMSSM specific contribution 5mh\rad arising 
because the PQ symmetry makes the lightest neutralino get a relatively small mass.^ This 
positive contribution is insensitive to tan/3, and becomes important for a small value of 
/ieff/AfsusY contrary to dmh\ 

mix- 



^See also for a similar discussion in a singlet extension of the MSSM having a relatively light 
neutralino. However, in our situation, the LEP bound on the invisible Z-boson decay width excludes large 
values of A at large tan /3. 
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Figure 1: The mass of the SM-Hke Higgs boson for Msusy = 1-5 TeV in the PQ-NMSSM. Here 
A is taken to be the maximum value satisfying the bound Tz-^j^o^^o < 2 MeV. The sohd curve in 
the left panel is the upper bound on rrih for iicS = 240 GeV (green) and fics — 420 GeV (blue), 
which is obtained for no mixing case Smhlmix = 0. The gray line is the MSSM value of m/j. In 
the left panel, the dashed lines are for Xt — while the solid ones for Xt = 6. Meanwhile, the 
right panel shows for Xt = 6, n^e = 420 GeV and a given value of Ax', the black curve is for 
A\ — 1.2MsusY while the red one for A\ ~ 2.4AfsusY- 



In Fig. I, we show the upper bound on irih in the PQ-NMSSM, which is obtained 
taking the maximum value of A allowed by the constraint ( ^.20| ) from the Z invisible decay. 
Here we have solved the RG equation to get Xh at the weak scale, and restricted A to be 
less than unity as would be necessary to maintain its perturbativity up to Mgut- In the 
left panel, the value of rrih is shown for 6mh\mix = 0. While the MSSM generates nih ~ 115 
GeV (128 GeV) at tan/3 > 10 for Xt = (6) and Msusy = 1-5 TeV, the additional 
contribution from A can lead to mh larger than 115 GeV also at low tan/3 as in other 
NMSSM models. The loops of stops involving At can further increase m/j. The maximum 
comes at Xt = 6. It is also important to note that, when ^gff ^ 400 GeV, rrih can be raised 
by a few GeV from the MSSM value even at large tan /3 owing to the PQ-NMSSM specific 
contribution 5m/i|rad- On the other hand, the right panel shows the value of nih for a given 
value of Ax. For Ax > 2^eff) the mixing effect is suppressed only at some limited region 
of tan /3. In the figure, we consider values of tan /3 giving (2/Xeg — Ax sin 2/3)^/m| less than 
0.1. 

3.2 Phenomenological aspects 

Since there appears a light neutralino as a consequence of the PQ symmetry, the Higgs 
boson h can invisibly decay into pairs of the lightest neutralino. This process is mediated 
by the Yukawa interaction /i-XiXii which is generated due to the mixing between S and 
^ and has a coupling given by 

..,0,0 = -V2X (<° sm/3 + <2 eos/3) . ^''jj"^^ (3.11) 

where the approximation is valid for small Xv / /igfj . The above coupling becomes negligible 
at large tan /3. If dominates, such non-standard invisible decay would make the Higgs 
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discovery at hadron colliders much more difficult. 

Recent LHC data have excluded the Higgs boson with SM properties in the mass range 
between 141 GeV and 476 GeV at the 95% confidence level ||2^. For h with mass lighter 
than 141 GeV, the main processes for its decay are h ^ bb and h — t- WW*, ZZ* pO| . 
The Higgs boson h in the PQ-NMSSM, which would have a small singlet component for 
/igjj/m| <C 1, can decay through a non-standard mode h — )• XiXi [01 ■ If it is kinematically 
accessible, the process h — )• XiXi takes place with the relative decay strength 

^W?x? ^ 1 (yj^y^ f AV f^^y /200GeVy 

^h^bi 3\m,/vJ \0.8j V 0.6 M /^cff y ' ^ ' 

r,^m 32vr2 /^/.xgx^y ... / 0.3 W A y /sin2/3y / 200GeV y 

r,^w* ~ 3i?(x) V <?w " "yRix)) \o.8j V 0.6 M /^eff y ^ 

where we have ignored the masses of the final states, and R{x) is defined by 

, 3(1 -8x + 20x2) /3x-l 
R[x) = —J- arccos ' 



(4x - 1)1/2 2x3/2 
-i-— ^(2- 13x + 47x2) - ^(1 -6x + 4x2)lnx, (3.14) 

with X = My^/mf^. The decay rate for the process h — )• ZZ* is similar to Fh^wiY*. When 
h has a sizable singlet component, the Higgs decay width for each process is modified, but 
the ratio between decay widths remains the same up to small correction arising due to 
that Uhxlxi receives contribution not only from H^H^ but also from SH^H^. The Higgs 
invisible decay to neutralinos would not dominate the SM decay processes either if A is 
strong enough to make x? heavier than m/j/2: 

^-°-^Hl30Giv) (200Giv) (ii^J ' 

or if A is small enough to suppress the Yukawa coupling of to the Higgs boson: Br(/i — 
XiXi) is less than 0.5 for 



1/2 



(3.16) 



Here we have naively estimated the value of A required for m^o > m/i/2 by taking the 
leading term in ( ^.14f) , which is expanded in powers of A^zj^/^^g. It is interesting to see 
that, in a low tan/3 region, the Higgs boson decays mainly through the invisible channel 
h — 5- XiXi for ^ 0.4 and /ieff ^ 400 GeV. A large /ieff can weaken this decay mode, 
but would lead to large mixing between and S. On the other hand, for tan j3 > 10, 



the constraint from the invisible Z-boson decay ( 2.20 ) requires A < 0.36 x (/icff/200GeV). 
Thus, in this case, Br(/i — )• x^/ti) cannot be larger than 0.5 if /igfj is smaller than about 
360 GeV. 

Fig. ^ shows the branching ratio for the invisible Higgs decay h — )■ XiXi iii the (tan /3, A) 
plane. In the figure, the yellow region is excluded by the experimental bound r^_^^o^o < 2 
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Figure 2: The branching ratio of non-standard mode h -> XiXi for ficS = 240 GeV (left) and 
HeS = 420 GeV (right). The dashed red hne is the contour for Br(/i XiXi)- We also show a 
contour plot for m^o larger than which is given in blue. The yellow region is excluded by 

the experimental bound on the invisible Z decay width. Meanwhile, the black contour shows the 
value of mil obtained for (5m/j|niix = in the case with Msusy = 1-5 TeV and Xt = 6. The Higgs 
boson has a mass larger than 115 GeV above the dashed black line. In the figure, masses are given 
in the GeV unit. 




1234512345 
tan/S tan/S 



Figure 3: The parameter region where Z — )■ xlxi is kinematically forbidden. The blue contour 
shows the value of m^o for /ies = 120 GeV. We also show the value of irih by a black contour for 
the case with Msusy = 1 TeV, Xt = 6 (left) and Xt =0 (right). The masses are given in the GeV 
unit. The dashed red line is a contour for Br(/i — > XiXi)- Above the solid red line, the lightest 
neutralino obtains a mass larger than mh/2. 



MeV, and the black contour is the value of nih obtained for Msusy = 1-5 TeV and Xt = Q 
with 5m/i|mix = 0. Here we have taken /iefr = 240 GeV for the left plot, and n^s = 420 
GeV for the right plot. Since there is an extra contribution to m/j from A as ( p.Sp , the 
Higgs mass can be raised above 115 GeV at low tan/3. Notice also that the contribution 
(^w-hlrad raises by a few GeV even at tan/3 > 10 compared to the MSSM value that 
corresponds to the A = case. The process h — t- XiX? can be the dominant mode of the 
Higgs decay at tan/3 < 10, but is suppressed at large tan/3. 
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On the other hand, for small /igff less than \v, the lightest neutralino is a sizable 
mixture of S and at low tan/3. In this case, it can acquire a mass larger than Mz/2 
so that the invisible decay Z — >• XiXi is kinematically forbidden. For some parameter 
region, it is possible for Xi to get a mass even larger than mh/2. Otherwise, h would decay 
dominantly through the invisible process h — )• XiXi low tan /3 because Xi has a sizable 
Higgsino component. One should also note that the contribution from A can raise m/j well 
above 115 GeV for MsusY ^ 1 TeV at tan/3 < 3. In Fig. ^, we show the region of A and 
tan/3 where Xi ticis a mass larger than Mz/2. The Higgs boson mass is also shown for the 
case with MsusY = 1 TeV. 



4. UV completion 

When one considers GUT models to UV complete the PQ-NMSSM, an important issue is 
how the GUT partners of the MSSM Higgs doublets, which also carry a PQ charge, acquire 
heavy masses. We point out that the PQ-NMSSM can emerge as a low energy effective 
theory of a missing-partner model for super symmetric SU(5) GUT [^].^ The missing- 
partner model has been considered to explain a large mass splitting of the SU(2) doublet 
and color triplet Higgses. The idea is to introduce higher dimensional representations 
that contain Higgs triplets but no doublets. Then, if a mass term H^H^ for 5 + 5 Higgs 
multiplets is absent in the superpotential, the doublet-triplet splitting can be achieved 
without fine-tuning from the interactions of and with the higher dimensional Higgs 
multiplets. This is possible because the superpotential is not renormalized in perturbation 
theory. 

The Higgs sector of the original model consists of the chiral multiplets, S(75), 9 (50) + 
^(50) and H(5) + H{5). To incorporate the PQ symmetry without spoiling the missing- 
partner mechanism, we modify the model by introducing three pairs of 50 + 50 and 5 + 5 
multiplets, and also three SU(5) singlets: 

S(75), 0,(50) + ^,(50), Hi{5)+Hi{5), Xi{l), (4.1) 

where i = (1,2,3), and C/(1)pq charges are assigned as 

S(0), Xi{q), X2{-3q), X^{2q), 

ei{-p)+h{p), e2{-3q-p) + e2{3q+p), 03{-q-p) + e3{q + p), 
Hi{p)+Hi{-q-p), H2{3q+p)+H2{-p), H^iq + p) + HsiSq - p). (4.2) 

This model seems similar to the minimal model with [/(1)pq considered in but it turns 
out that more than two pairs of chiral multiplets are needed to obtain the PQ-NMSSM as 



®A 5 dimensional SU(5) unified theory can also yield the PQ-NMSSM when compactified on /{Z2 x Z2) 
orbifold |3^] . For instance, one can introduce a pair of Higgs hypermultiplet H + % which form 5-1-5 
representation of SU(5) and carry a PQ charge —1. Then, doublet-triplet splitting is achieved taking the 
orbifold projection such that H (H) has a SU(2) doublet Higgs chiral multiplet transforming as (+,+) 



under Z2 x Z2 and a triplet Higgs with (+, — ). This also leads to the terms (2.1) for the PQ-NMSSM below 
the compactification scale. 
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an effective theory below Mqut- The missing-partner mechanism is implemented by the 
following PQ-invariant superpotential terms 



W = ^MT¥(e2) + ^aVi{Y?) + hA^Hi + cA^Hi+i + MA0i, (4.3) 

with the identification ^4 = Hi. As in the original model, the vacuum expectation value 
of T, breaks SU(5) to the SM gauge groups, and gives rise to the triplet mass terms. This 
becomes clear after integrating out the heavy triplets in 6i + 6i, which leads to 

W,s = MfH^H^_,,, (4.4) 

where M? ~ M^^r^/ Mi, and denotes the color-triplet from Hi. Because the 50 repre- 
sentation does not contain Higgs doublets, no doublet mass terms are generated from the 
superpotential ( [4. 3D , and the three pairs of doublet Higgses remain massless. Mass terms 
for these doublet Higgses arise from 

W = XiX.HiH, + -^XfX2, (4.5) 
Mpi 

where we have chosen a basis of X3 such that XfH^H^ is removed in the superpotential. 
Including soft SUSY breaking terms for the gauge singlet scalars 

^2 I |2 I / /I ^ " 



- £soft = m'xjXil' + \^A^^XtX2 + h.c. j , (4.6) 

the second term in the above superpotential fixes Xi 2 at 

Mp,(-m^jV^ MpiMsvsY 
\Xi^2\ w ^ '■ ^ , (4.7) 

thereby leading to that J7(1)pq is spontaneously broken around 10^^ GeV for ^ ~ 1, and 
there appears the axion which is a mixture of arg(Xi^2)- Here we have assumed m'j^_^ < 0. 
Hence, the doublet Higgses in -ffi,2 + -ffi,2 obtain large masses from the vacuum expectation 
value of Xi^2, respectively. On the other hand, the doublet Higg ses in ^3 -t- H^ remain 
massless until X3 acquires a nonzero vacuum expectation value. 



Notice that the PQ charge assignment (4^) allows direct mass terms HiHi-^i and the 



Yukawa term X1X2X3 in the renormalizable superpotential. These terms should be absent 
in order for the missing-partner mechanism to work and for only one pair of Higgs doublets 
to remain light. Once we do not put these superpotential terms, radiative corrections will 
not change the situation owing to supersymmetry.^ It is also important to note that the 
model possesses two global U{1) symmetries associated with the independent charges p 
and q. To eliminate one of them, as was considered in [^, we introduce three right-handed 
neutrino multiplets A^(l) that implement the conventional see-saw mechanism through the 



^One can assign a different PQ charge to X2 and the Higgs multiplets to forbid a superpotential term 
X1X2XS. Then, other mechanism is needed to fix the PQ breaking scale because the term XfX2 in the 



scalar potential (4.6) is not allowed 
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superpotential terms NLH^ + XiNN with i = 1 or 2. The PQ charges are then fixed 
as 5p = —8q when the Majorana masses for arise from XiNN, and bp = —12q if one 
instead chooses X2NN. 

It now becomes apparent that the missing-partner model with the superpotential terms 
( |4^ ) and (|]|) leads to the PQ-NMSSM. The doublet Higgses in H3 + H3 correspond to 
the ordinary MSSM Higgses, while Xi^2 and X3 play the role of X and S, respectively. In 
the model, the higher dimensional operators 



where ttiq ~ ^SUSY '^i.s ~ 1 and ^ ~ 1. Then, X3 is naturally expected to get a vacuum 
expectation value around the weak scale for MsusY ~ 1 TeV. However, if PQ-breaking 
mass terms OiOj with i ^ j are present, the loops of heavy triplets would generate large 
tadpoles for X^. This implies that f7(l)pQ is crucial to avoid the tadpole problem. The 
PQ symmetry plays an important role also in suppressing dangerous higher dimensional 
operators leading to too rapid proton decays. The triplet Higgses mediate dimension 5 
operators violating the baryon number [^3|, which carry a nonzero PQ charge and therefore 
are further suppressed by a small factor XiX2/(M|^M|) ~ Msusy-^p«/(-^i-^2) compared 
to those in the minimal SU(5) GUT model. 

Let us finally discuss the difference from the model of That model contains two 
pairs of Higgs doublets Hf + Hj and Hj + H'j: which are vector-like also under C/(1)pq. 
One pair of them becomes heavy through PHfH'j for P being a C/(1)pq breaking gauge 
singlet field, and the other remains light. If one introduces an additional singlet P' having 
a term P'HjHf in the superpotential, P' necessarily carries a PQ charge such that PP' is 
invariant under f7(l)pQ transformations. Thus, even if one omits PP' in the superpotential, 
a Kahler potential term PP' would induce too large tadpole term for P'. This makes it 
difficult for P' to play the role of 5 in the PQ-NMSSM. 

5. Conclusions 

Extended to incorporate the PQ mechanism solving the strong CP problem, the NMSSM 
becomes compatible with the grand unification since the PQ symmetry forbids large tad- 
poles for the SM singlet S to be generated from loops of heavy fields coupling to S. Another 
important property of the PQ-NMSSM is that all the mass parameters are determined by 
the SUSY breaking scale MgusY and F^/Mpi with Fa being the axion decay constant. 
Thus, the electroweak symmetry breaking is naturally achieved at the correct scale. Fur- 
thermore, the model can avoid the domain wall problem in the presence of the PQ mes- 
sengers. 

An important consequence of the PQ symmetry is that the lightest neutralino is 
singlino-like with a small Higgsino admixture, and is relatively light compared to other 




(4.8) 



can generate an effective tadpole term for X3 as 



(4.9) 
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sparticles. The Higgsino component is determined by the couphng of S to the Higgs dou- 
blets, which is constrained by the LEP bound on the invisible Z-boson decay width. This 
constraint becomes severe at large tan /3. Meanwhile, the SM-like Higgs boson decays 
mainly through the conventional decay modes at large tan /3 and in a portion of parameter 
space for small values of tan /3. The decay of the Higgs boson into a pair of the lightest 
neutralino can be the main mode at low tan f3, for which case the Higgs search at colliders 
will be modified. Also important is that the SM-like Higgs mass receives an additional 
positive contribution from the loops involving the singlino Yukawa coupling. This PQ- 
NMSSM specific contribution can lead to a significant increase of the Higgs boson mass by 
a few GeV even at large tan /3 compared to the MSSM. 

We found that the PQ-NMSSM is realized as a low energy effective theory of a missing- 
partner model for supersymmetric SU(5) GUT with the PQ symmetry, which solves the 
doublet-triplet splitting problem and the proton decay problem. It is interesting to note 
that such a UV completion achieves the relation Fa ~ \/MsusY-^Pi- Hence, all the mass 
parameters of the resulting PQ-NMSSM have values of the order of MsusY- 

Note added 

After submitting the manuscript, the ATLAS and CMS collaborations at the LHC 
reported their updated results in the Higgs search which may indicate a SM-like 
Higgs boson with mass around 125 GeV. To explain a 125 GeV Higgs mass within the 
MSSM, we need large stop mixing or heavy stops with mass larger than about 10 TeV. 
The PQ-NMSSM improves the situation because the Higgs mass receives an additional 
positive contribution, which can be of a few GeV even at large tan/3. 
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A. Mass matrix 

In this appendix, we present the tree-level mass matrices for the neutral scalar fields. After 
rotating the upper left 2x2 submatrix of the mass matrix for the CP even scalars, one 
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obtains 

(M|)ii = M| cos^ 2/3 + A^w^ sin^ 2/3, 

(M^)22 = ^ + (M|-AV)sm2 2/3, 
sm 2/5 

(M^)33 = m| + AV, 
(M^)i2,2i = ^(M|-AV)sm4/3, 
(M^)i3,3i = A?;(2/ieff-^Asm2/3), 

(M^)23,32 = A^;AACOs2/3. (A.l) 
The mass matrix for the pseudoscalar fields is given by 

^ 6eff cot /3 6eff Aw^A COS /3 \ 

6cff 6efftan/3 Au^dAsin/? . (A. 2) 

y \vA\ cos /3 Af sin /3 m| + A^w^ y 

It is easy to see that there are one massless mode, which is absorbed into gauge boson, and 
two massive CP odd scalars: 



Using the stationary condition (p.7p, one can find that = if k = 0. 



B. Global structure of the Higgs potential 

In the PQ-NMSSM, where the Higgs sector is extended to include the singlet S", the Higgs 
potential may develop another minimum away from the weak scale. The model parameters 
are constrained to avoid such a minimum since it would generally appear at a field value 
similar to or larger than MsusY and thus be deeper than the electroweak vacuum. For the 
case with m| > and A < 1, we shall show that the region A^(|//°p + |-f^°P) ~ ^SUSY 
somewhat near the D-flat direction is potentially dangerous for much of the parameter 
space. This argues that the condition to avoid a deeper minimum can approximately be 
examined by looking at the shape of the potential along the D-flat direction. 

After integrating out S by the minimization condition, the Higgs potential reads 



2A202 + ^2 

+ 2 {rn\^^ sin^ Q + m^^ cos^ Q - Xm^ sin 26') + constant, (B.l) 

where = ^(\)smQ and = cos d with < < 7r/2 and < (/>. At very large 
values of the first term becomes dominant and lifts the potential along the 0-direction. 
It is also straightforward to see that = when = or when 

sin2 20 + cos2 20 j {2\^<^'' + - A:m| {2\^^^ + m|)' 

+ 2 r^^^sin20- AS^mg^j =0, (B.2) 
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where is a function of 9, 

( A? \ n + q'"^ o 'm'jr sin^ 6 + mi, cos^ 6 - Xml sin 26 

k= 1 + -1^ sin^ 26+ ^ cos^ 29 - 2—^ ^ • (^-3) 

\ 2m^ J rrig 

The above relation shows that the potential can have at most one local minimum at (p ^ 
along the c/i-direction for a given 9, which would appear at 2}?(f)^ + m| ~ -^SUSY- 
other hand, along the angular direction, the slope of the potential vanishes when 

/ (7,2 

sin 26* = - (^r^ sin 261 + COS 261, (B.4) 



for which 



Here r is defined by 



dlV = 4 cos 2^ + ^ sin 29^ Af (/>2 tan^ 29. (B.5) 



and 2 ^'^^ given by 



A 2 2 2 

both of which are generally of ©(Mgygy), and positive. For r > 0, one can find (i) = 
can have a solution at tan0 > 1 with a positive curvature d'^V > 0, implying that there is 
only one minimum along the angular direction for a given (p, and (ii) for large values of (p, 
0^ ^ A^ ~ -^SUSY' ^ rninimum along the angular direction is located near tan0 = 1, i.e. 
near the D-flat direction. 

Let us examine further the case with r > and m|, ~ -^SUSY- Note that r is positive 
at 2X'^(j)'^ + m| ~ ^susY if ^ ^ 0-5 for small Ax, and if A < 1 for Ax ~ MsusY and 

~ AfsusY- At the electroweak vacuum, which lies at cp'^ ~ <^ -^SUSY ^^'^ ^ = /3, 
the condition ( |B.4| ) gives 

tan 2/3. (B.8) 

^1 



The extremum condition {B.4) can then be written 



^_r^ tan^g ^ ) tan 20 tan 2/3, (B.9) 



Af tan2 6' + l, 

for (/>^ < ^SUSY- This tells that a minimum along the angular direction arises at 1 < 
tan0 < tan/3 for a given (p, and approaches the D-flat direction, tan0 = 1, as 4> increases. 
Thus it is useful to first analyze the potential along the D-flat direction though the actual 
another minimum, if exists, appears somewhat away from the D-flat direction. 

On the other hand, in the case with r < 0, the potential is minimized along the angular 
direction at tan^ > tan/3 or at tan 9 < 1 for </>^ > -^SUSY- former case, making the 

potential develop no other minimum in the region near the iD-flat direction is not enough 
to guarantee the absence of a deeper minimum. 
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C. RG running equations 



For the low energy effective theory ( |3.1[ ) below MsusY, the RG running equations for the 
Yukawa couplings read 

8vr fJ--r- = + 4y„ + -y^ - —g^ - -g^ , (C.l) 



at the one-loop, and those for SM gauge couplings are 



2 d 1 103 1 , 17,, , 2 

dfi gf 30 10 30 5 



dfi g2 o 6 



87r^ = 8-9{fi-mt), (C.2) 



with (7i = y 5/3(7' and 32 = 9- Here 9{x) = 1 for x > and 9{x) = for x < 0. 
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